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Rapid Modeling Based on Parametric Interface Representation 
-Final Report for NAG-1-1670 - 

Background 

This report summarizes the work of the graduate student, Mr. Ollie J. Rose, for the period 
from January 1, 1996 through August 1, 1997. Mr. Rose began this research activity under Dr. Mo- 
hammad Aminpour. However, Dr. Aminpour left ODU at the end of the Fall 1995 semester. At that 
time. Dr. Norman F. Knight, Jr. agreed to direct Mr. Rose’s research work and Dr. Osama A. Kandil 
agreed to represent the University as the PI on the grant. The grant remained open under a no-cost 
extension under August 1997. There were approximately $ 1 1 ,000 remaining in the grant for student 
support which was paid to Mr. Rose as a graduate research assistantship stipend along with additional 
Departmental funds. Mr. Rose is nearing completion of his dissertation which is anticipated by De- 
cember 1998. At that time, a copy of the dissertation will be forwarded to the contracting officers 
technical representative. Dr. W. Jefferson Stroud, of the NASA Langley Research Center. 

Introduction 

Over the past few years, the development of the interface technology (see Ref. 1-1 1) has pro- 
vided an analysis framework for embedding detailed finite element models within finite element 
models which are less refined. This development has enabled the use of cascading substructure do- 
mains without the constraint of coincident nodes along substructure boundaries. The approach used 
for the interface element is based on an alternate variational principle often used in deriving hybrid 
finite elements. The resulting system of equations exhibits a high degree of sparsity but gives rise 
to a non-positive definite system which causes difficulties with many of the equation solvers in gen- 
eral-purpose finite element codes. Hence the global system of equations is generally solved using 
a decomposition procedure with pivoting. 

The research reported to-date for the interface element includes the one-dimensional line 
interface element (e.g. Ref. 7) and two-dimensional surface interface element (Ref. 9). Several 
large-scale simulations (e.g. Ref. 4 and 11), including geometrically nonlinear problems (Ref. 10 
and 11), have been reported using the one-dimensional interface element technology; however, 
only limited applications are available for the surface interface element. In the applications reported 
to-date, the geometry of the interfaced domains exactly match each other even though the spatial 
discretization within each domain may be different. As such, the spatial modeling of each domain, 
the interface elements and the assembled system is still laborious. 

The present research is focussed on developing a rapid modeling procedure based on a para- 
metric interface representation of independently defined subdomains which are also independently 
discretized. The geometric definition of the interface element, either one-dimensional or two-di- 



mensional, is automatically made using a parametric representation requiring very little input from 
the analyst. In many cases where detailed finite element models are desired locally within a given 
structure, this rapid modeling approach can expedite the modeling task and lead to faster turn around 
between design and analysis. 

The dissertation will present and describe this rapid modeling approach. The approach is 
an integral part of the interface technology and provides a new capability not described or available 
in the existing interface element work. Applications of a new, spatially curved (non-planar) inter- 
face and a general two-dimensional surface interface will be discussed and demonstrated. 

Approach 

The approach used to define the interface of two or more independently defined domains 
includes an assessment of the domains to be interfaced and the specification of weighting functions, 
if desired, to guide the creation of the interface element geometry in a desired manner. Node defini- 
tion for the interface can include nodes from any of the adjacent subdomains or alternatively may 
be specified as an independent curve or surface. Consider, for example, a general spatially curved 
interface element. Here the points along the edges of the domains to be connected are determined 
automatically based on connectivity information. Using this information, the geometry of the inter- 
face element is created by using a least-squares combination of cubic spline basis functions with 
appropriate weighting functions, if specified. In the case of a space curve, coordinates on the inter- 
face element are represented parametrically in three dimensions as follows: 
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In this representation, Bj(s) are the single-parameter cubic spline basis functions, s is a parametric 
coordinate along the one-dimensional interface, the a x (i), a y (i), and a z (i) are coefficients deter- 
mined by the least-squares process and N is the number of basis functions (N depends on the knot 
distribution and order of the spline polynomials). For the case of two-dimensional surface inter- 
faces, the interface coordinates are represented in a similar manner, but the basis functions become 
tensor products of linearly-independent, single-parameter basis functions to form a two-parameter 
representation. 

Preliminary Results 

Results are presented for the one-dimensional interface element based on the proposed rapid 
modeling procedure. The features and capabilities of this rapid modeling procedure are demon- 



strated using two-dimensional and three-dimensional domains which are independently defined 
and discretized. The two-dimensional problem shown in Figure 1 illustrates the common problem 
of a circular cutout in a rectangular domain. Previously, the interface element would have to be de- 
fined at the outset by the analyst. Using the present approach, one subdomain could be an annular 
region defined by two concentric circular arcs, while the other domain is a rectangular domain with 
a square cutout, as shown in Figure 2. For the case of uniform thickness and material properties, 
the rapid modeling approach then creates an interface boundary and automatically redefines the sub- 
domain discretizations in such a way as to eliminate any gaps that may have existed in the original 
subdomain geometries. An example of independent grids for the plate with cutout and the resulting 
model after the remeshing step are shown in Figure 3. 

The robustness of the parametric interface is dearly evident when applied to curved domains 
such as a cylinder with a circular cutout. An example of this is given in Figure 4 which shows the 
independent domains and adapted spatially curved interface, indicating the advantages of the present 
method. 

The new one-dimensional parametric arc interface representation is being extended to two- 
dimensional parametric surface representations. This extension is nearly developed and represents 
a significant advancement in interface element technology. An example of independent grids for 
the plate with cutout and the resulting 3D model after the remeshing step are shown in Figure 5. 

Summary 

The dissertation will describe the parametric interface representation process and the least- 
squares, cubic spline functions used. A brief description of the interface element technology will 
be included for completeness. Applications of this rapid modeling approach will be presented for 
flat two-dimensional plane stress problems and three-dimensional elasticity analyses to illustrate 
the robustness of the present approach for modeling and simulation. Displacement and stress results 
will be given for representative examples cases. 
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Figure 1: Circular Cutout in a Rectangular Plate 



Figure 2: Independently Modeled Subdomains 



Independently Modeled Domains 
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Figure 3: Automatic Interface Modeling of a Rectangular Plate with Circular Cutout 



Figure 4: Automatic Interface Modeling of a Cylindrical Panel with Circular Cutout 
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Status Report for NASA Grant NAG1-1670 
I. Work Completed \2J/5/^S 

The proposal for Grant NAG 1-1 670 was written for a three year period (please 
see the Duration Section in the proposal, a copy of which is attached). However, 
the grant was budgeted only for the first year. Therefore, the work reported here 
pertains only to this first year and not the entire work stated in the proposal. In 
this first year phase of the grant the following has been completed. 

A. In the course of studying and testing the 1-D interface element (as a 
prelude to the 2-D work) errors were encountered in the code for built- 
up structures and curved interfaces. These errors were corrected and 
the benchmark Boeing composite crown panel was analyzed 
successfully. 

B. Studied the NURBS (Non-Uniform Rational B-Splines), cubic splines, 
and cubic B-splines and concluded that the cubic B-splines are the 
state of the art and best suite our application. 

C. When analyzing the Boeing crown panel, it was observed that 
coordinate data points representing the interface contained small 
random error (noise) which in turn caused some noise in the results. 
This noisy coordinate data will also be present in the 2-D interface, 
perhaps causing more severe error in the results. To construct a 
smooth function from the noisy data a least squares approach 
combined with cubic B-splines for the planar curve in the 1-D case was 
developed and tested with success. The work on planar curves should 
be generalizable if need for space curves arises in the future. The 
theoretical formulation is designed such that the 2-D surface 
representation would be realized as a tensor product of the 1-D work, 
thus capitalizing on the 1-D work. 

D. Following the successful formulation of a discontinuous 1-D beam 
element, preliminary investigations for the formulation of discontinuous 
2-D shell and 3-D solid elements were conducted. It turns out that the 
2-D hybrid shell elements can be reformulated (some programming is 
required) to account for discontinuity across the shell element surface 
(applicable in the cross-surface method mentioned in the proposal). 
However, for discontinuous loading on the shell element edges and on 
the 3-D solid surfaces, although a formulation is workable using 
classical elasticity solutions, the formulation is not as straight forward 



and the accuracy of the results is not as clear. More investigation is 
necessary for these cases. 

Some key features of the technology developed so far are summarized below. 

1 . The basis functions are piecewise cubic polynomials defined over the 
entire domain and are continuous through the second derivative. 

2. The basis functions are of compact support, i.e., they are non-zero over 
only a 5-knot consecutive subset of knots. This implies that solution 
matrices derived for fitting data are pentadiagonal. 

3. The knot set and corresponding basis functions are defined 
independently of the data ordinates. There is an algorithm for automatic 
knot distribution definition which provides nearly equal number density 
of data points in the knot partition domain. 

4. The data can be given in local or global coordinates, but the fitting takes 
place locally for more accurate data representation. 

A sketch of a typical 1-D interface between incompatible grids is given in diagram 
A. Examples of the 1-D "noisy data" representation capability are shown in 
Figures 1-6. Figure 1 shows basis functions for a non-uniform knot set of nine 
knots partitioning the domain [0,1]. Basis function number four is shown alone for 
clarity in Figure 2 plotted over the entire domain. The region of compact support 
is clearly seen. In Figure 3 is plotted a circular arc section of 10% maximum 
height together with the least-squares fit using the previously shown 9-knot basis 
functions. The vertical scale is exaggerated by a factor of about 10, and the 
excellent quality of representation can be seen. The same circular arc data with 
considerable random error corrupting the points is shown in Figure 4, together 
with the least-squares spline fit. Again, the scale is exaggerated and the basis 
splines represent the data nicely. 

Figures 5 and 6 show circular arc augments of 90 degrees and 135 degrees in a 
global coordinate system. The data are fitted locally with the knot set spanning 
the chords. This technique of local fitting generally assures a single-valued 
function and improved accuracy. 



II. Work in Progress and Planned 

A. Two Dimensional Interface Representation 

1 . An algorithm is being developed to select a local coordinate 
system to serve as a 2-D parameterization space and projection 
plane. Various issues and options for this algorithm are 
currently being studied and assessed. 

2. The 1-D knot partitioning and distribution algorithm is being 
generalized for 2-D partitioning and knot selection in the local 
projection plane. The algorithm will be such that nearly equal 
number density of projected points occur in each knot-partition 
cell. 

3. Development of 2-D basis functions as tensor products if 1-D 
basis functions. 

4. Algorithm for least-squares representation of the interface using 
2-D basis functions. 

B. Two-Dimensional Interface Element 

1. Formulate equations coupling displacements and tractions using 
a hybrid variational method similar to the 1-D formulation. 

2. Generalized stiffness matrices resulting from the above 
mentioned formulation contain sub-matrices referred to as the M 
and G matrices. These submatrices are obtained by integration, 
over the interface, of quantities related to traction and 
displacement. These integrals are to be computed numerically 
using the bi-cubic B-splines as interpolation functions for 
integration. 

C. The formulation and implementation of consistent load calculations for 
the 3-D hybrid solid elements in COMET-AR. These load calculations 
are necessary to enable the use of the hybrid solid elements in 
COMET-AR in conjunction with the 2-D interface element. 



D. Benchmarking and Demonstration Applications 


1 . The 2-D interface technology is expected to be verified by solving 
a beam with a 2-D planar interface and comparing to existing 
solutions. 

2. A second problem of suitable complexity and having a curved 
interface is to be solved. The exact features of the proposed 
problem have not been decided upon, but an example of the type 
of problem envisioned is a ball and socket for two or more 
interfaces resulting from relative movement between the ball and 
socket. The envisioned methodology will be restricted to 
movement increments at least one grid line in one of the grid 
directions. A pictorial sketch of an expected typical 2-D interface 
is shown in diagram B. 



Diagram A: Curved 1-D Interface 






















REPRESENTATION 
RADIUS 5 




Duagram B: Curved 2-D Interface 



A Two Dimensional Interface Element for Coupling of Independently 
Modeled Three Dimensional Finite Element Meshes and Extensions 
to Dynamic and Non-Linear Regimes 


Background 

A newly developed one dimensional interface element 1-4 was shown to be a 
successful method for coupling of independently modeled two dimensional (e.g., 
plate and shell) finite element meshes. The use of the interface element eliminates 
transition modeling and greatly simplifies the modeling efforts. As such this 
interface element may be used for effective modeling of aerospace structures by 
connecting independently modeled substructures as well as independently modeled 
global and local subdomains. 

This interface element was based on the hybrid variational principles and was shown 
to be able to recover the constant and linear states of strain and curvature with no 
error (provided that the finite elements used in the analysis are capable of doing 
the same). The exact recovery of these states of strain and curvature is essential 
if one hopes to converge to the correct solution as the associated finite element 
meshes and the interface elements are refined. Subsequently these tests may be 
considered as patch tests for the interface element. The interface element based on 
the hybrid variational principle possesses the following four inherent properties. (1) 
The balance of energy is preserved across the interface (i.e. the constraint integrals 
introduced into the formulation do not contribute to the energy of the system). 
(2) Equilibrium is maintained across the interface (i.e., the sum of the forces at 
the interface is zero). (3) Displacement compatibility is preserved, in a variational 
sense, across the interface. (4) The forces at the interface are redistributed among 
various substructures in a variational sense that is consistent with the distribution 
of forces over the finite elements at the interface from each substructure. 

The first two properties are satisfied by any method that produces a correct result 
on a global scale. The last two properties, however, are not satisfied by most 
methods and as a result a high level of noise is usually introduced into the solution 
details (such as stresses) particularly near the interface region. In these methods, 
the property (3) above, is satisfied only at discrete points and as a result the forces 
at the interface are redistributed among various substructures in a manner that is 
not consistent with the distribution of forces over the finite elements at the interface 
from each substructure. These forces, however, are in an overall equilibrium 
(i.e., property (2) above) but do not have the “correct” distribution. The hybrid 
variational formulation, on the other hand, provides the “correct” distribution of 
forces along the edge of each substructure at the interface. This is the key to the 
success of the hybrid variational formulation. 



Also, the 1-D hybrid variational interface element formulation has recently been 
extended to develop an enhanced interface element that possesses cross surface 
capabilities 5 . While the original interface element is placed along substructure 
edges, this new enhanced interface element may be placed over a surface crossing the 
element sides. While the original element simplifies modeling and avoids transition 
modeling in substructuring and global/local analysis, this latter enhancement 
simplifies modeling and enables the user to perform component substructuring. 
That is, the components may first be modeled independently and then assembled 
together and not necessarily along their edges. However, this latter element, 
although successful, needs some more work before becoming fully operational. 

The interface element may also be utilized to perform “substructure adaptive 
refinement”. The major difficulty with the current adaptive refinement strategies is 
that they create distorted meshes which, in turn, deteriorate the performance of the 
underlying finite elements and create the so called mesh locking phenomenon. With 
the use of interface elements, however, it is natural to perform uniform refinement in 
each subdomain such that the underlying finite elements are not distorted beyond 
their original shapes and, hence, the performance of the finite elements are not 
deteriorated. This strategy may be referred to as a “adaptive-uniform refinement 
technique.” 

The major drawback of the interface element technology is that the final system 
matrix, although nonsingular, is an indefinite matrix. As such, a unique solution 
exists for the problem at hand; however, conventional positive definite fast solvers 
cannot be used for this method, and regular Gaussian elimination solvers are too 
slow and require a large amount of memory. There are two tentative remedies 
to this difficulty. 1) To develop a special sparse solver with full pivoting (single 
and double rows and columns) with out-of-core capability such as the one in the 
Boeing extended mathematical library and 2) to partition the system matrix and 
manipulate it to create smaller size positive definite matrices that can be solved 
using the available fast solvers. 

Proposed Research 

Since the one dimensional interface element has been a very successful method 
for coupling of shell structures, it is proposed herein to explore the possibilities of 
coupling of three dimensional finite element subdomains using a two dimensional 
interface element. This two dimensional interface element does not exist and its 
development requires research. Some type of surface spline has to be used for 
the development of this element in order to accurately match the surfaces of the 
independently modeled 3-D solids. Surface splines in conjunction with global/local 
analysis have already been used 6,7 and the technology is available to us with 
little effort. However, the performance of these surface splines in this kind of 
application has to be assessed and most likely the NURBS (Non-Uniform Rational 



B-Splines) surface design methods will have to be formulated and programmed in the 
development of this interface element. However, this requires some research and it is 
not as straight forward as the surface splines that have been used in conjunction with 
the global/local analyses. The Lagrange multipliers interpolation used in the 1-D 
interface element also has to be extended to the 2-D version. This also requires some 
research regarding the type and order of approximations used to ensure numerical 
stability and rank sufficiency of the final system matrix. This element would be 
very useful for 3-D/3-D and 2-D/3-D substructuring and global/local analysis in 
conducting local detail stress analysis of primary structural components. 

Also, since the enhanced version of the 1-D interface element for component 
substructuring has been successful, it is proposed herein that this concept to be 
also extended to 3-D structural applications which makes modeling of complicated 
structures very simple by independent 3-D component substructuring. 

It is also proposed that the extension of the 1-D interface element to dynamic and 
nonlinear regimes be studied. These seem to be natural extensions for the successful 

1- D interface element. Furthermore, it is proposed to study the extension of the 

2- D interface element to dynamic and nonlinear regimes as well. 

Finally, the issues discussed earlier in the background section regarding an efficient 
equation solver will be studied and the two possibilities discussed there will be 
explored to devise a satisfactory and workable solution method. 

Once the formulation and implementation is complete, several test cases will be 
studied to ensure accuracy and correctness of the formulation. These test cases will 
include various boundary conditions and loadings. 

Personnel and Resources 

Dr. Aminpour will be the principal investigator. He will work with two of his 
graduate students to accomplish this research. NASA computing facilities will be 
used and some time will be spent at NASA LaRC during the period of performance. 

Duration 

The above research is proposed to be performed over a period of three years (October 
1, 1994 through September 30, 1997) by the principal investigator on a part time 
basis during the academic year and 100% of his time during the su mm er and two 
, graduate students at 50% of their time throughout the period of performance. 
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Consistent Load Calculations for 3-D Hybrid Solid Elements in 
COMET- AR and Investigation of 2-D and 3-D Discontinuous 
Finite Elements to Be Used with Interface Elements 
Supplemental Task No. 1 for NAGl-1070 


Supplemental Task 1 

Supplemental task 1 will focus on the following two topics. 

(i) The formulation and implementation of consistent load calculations for the three-dimensional hybrid solid 
elements in COMET-AR. These load calculations are necessary to enable the use of the solid elements in 
conjunction with the 2-D interface element under development. This 2-D interface element, like it's 1- 
D counterpart is based on the hybrid variational principle 1-3 and is used for connecting independently 
modeled three-dimensional substructures. 

(ii) When connecting independently modeled substructures) the elements on the boundary of one substructure 
may cut across the faces of the elements on the boundary of other substructures creating a stress discontinuity 
in those elements. Current finite elements are not capable of dealing with these internal discontinuities. In 
references 4 and 5, two different methods were devised to deal with this situation for the 1-D interface 
element. Both methods resulted in satisfactory results. However, these methods are somewhat awkward 
to implement and the extension of these methods to the 2-D interface element is not very promising. It 
is proposed that a preliminary investigation into the formulation of discontinuous 2-D shell and 3-D solid 
elements be conducted. A formulation of a discontinuous 1-D beam element by the Principal Investigator 
has already been developed with very satisfactory results. 

Personnel and Resources 

Dr. Aminpour will be the Principal Investigator. He will work with one of his graduate students to accomplish 
this research. NASA computing facilities will be used and some time will be spent at NASA LaRC during 
the period of performance. 

Duration 

The above research is proposed to be performed over the period September 1, 1995 through January 11, 
1996 by the principal investigator on a part time basis and one graduate student at 50% of his/her time 
throughout the period of performance. 
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